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Abstract 

We consider critical dense polymers £(1,2). We obtain for this model the eigenvalues of the 
local integrals of motion of the underlying Conformal Field Theory by means of Thermo- 
dynamic Bethe Ansatz. We give a detailed description of the relation between this model 
and Symplectic Fermions including some examples of the indecomposable structure of the 
transfer matrix in the continuum limit. Integrals of motion are defined directly on the lattice 
in terms of the Temperley Lieb Algebra and their eigenvalues are obtained and expressed as 
an infinite sum of the eigenvalues of the continuum integrals of motion. An elegant decom- 
position of the transfer matrix in terms of a finite number of lattice integrals of motion is 
obtained thus providing a reason for their introduction. 

1 Introduction 

It is well established that a lattice approach to logarithmic minimal models Cp^q [4] can be 
realized in terms of indecomposable representations of the Temperley Lieb Algebra [3], in 
particular the integrability of these lattice realizations of logarithmic CFTs is proved by the 
existence of commuting families of double row A— tangles, the parameter of such a family 
being called the spectral parameter. 

In contrast with unitary minimal models, which are realized on the lattice for example by 
the RSOS models [19], the transfer matrix may exhibit a Jordan indecomposable structure 
for some choice of Cardy-type boundary conditions. Furthermore the Cp^q are not known so 
far to posses an elliptic deformation as the critical RSOS models. Defining for these models 
a scaling limit procedure leads to additional divergences due to the model being critical and 
therefore not being possible to define a limit in which a combination of the elliptic modular 
parameter and the system size generate a flow parameter driving an integrable perturbation 
of CFT away from the UV fixed point. Such a flow parameter would then act as a cut off 
regulating the divergences in the Thermodynamic Bethe Ansatz (TBA) integral equations. 
The perturbed CFT associated with the continuum limit of the elliptic deformation of the 
critical A3 model is for example the thermal 0i,3 perturbation of the minimal CFT A^3,4. 
In this paper we want to discuss in detail the model £1^2 (also called Critical Dense Poly- 
mers [5]) in view of the well known works of Bazanov, Lukyanov and A.B. Zamolodchikov 
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(BLZ)[1], and describe exhaustively how the local involutive BLZ charges of Conformal Field 
Theory arise in a number of expansions directly on the lattice. In particular we shall learn 
something along the way of the relation between critical dense polymers and symplectic 
fermions. 

It is well known [5] that the CFT corresponding to critical dense polymers has central charge 
c = —2. Such a conformal field theory is known to be logarithmic, these theories, in contrast 
with rational CFTs, can be realized by different models for the same value of the central 
charge and conformal weights. For example Hamiltonian walks on a Manhattan lattice [7, 8] , 
the rational triplet theory[9, 10, 11], symplectic fermions [12, 13], the Abelian sandpile model 
[14], dimers [15], the traveling salesman problem [16], branching polymers [17] and spanning 
webs [18] all share the same value of the central charge, which is —2. 

The layout of the paper proceeds by reviewing some common lore about the CFT corre- 
ponding to critical dense polymers. In section 2 the lattice model is introduced, the transfer 
matrix is explicitly built from the Boltzmann weights and the inversion identities and selec- 
tion rules are also discussed. 

In section 3 we derive the TEA equations for the model and deal with its continuum limit. 
The eigenvalues of the BLZ involutive charges are obtained by expanding the eigenvalues 
of the continuum scaled transfer matrix. In section 3.3 a new result is obtained, that is 
after having identified the involutive charges we are able to perform a 1/A^ expansion for 
the eigenvalues in which the conserved charges explicitly appear. Such an expansion is then 
manipulated to obtain an alternative form that provides the eigenvalues for the lattice invo- 
lutive charges. In this new framework the eigenvalues of the transfer matrix are expressed 
in terms of Bell polynomials, and the inversion identity itself is expressed in terms of these 
polynomials. 

In section 4 these results are extended to the transfer matrix itself and the iV— tangles cor- 
responding to the lattice involutive charges are explicitly built in terms of the Temperley 
Lieb algebra, thus providing a reason for the long calculation of section 3.3. 
In section 5 we describe the relation of the model with symplectic fermions. We give a 
description of selection rules for (r, s) boundary conditions which is analogous to the lattice 
one, we decompose all the characters in terms of characters of certain fermionic modules 
built over the Virasoro algebra. And finally we describe the Jordan decomposition of the 
continuum transfer matrix corresponding to modules with the same conformal weight but 
different (r, s). It is shown how (1, s) modules correspond to diagonalizable transfer matri- 
ces in agreement with the lattice behavior of the model, on the other hand for r 7^ 1 it can 
happen that the transfer matrix exhibits a nontrivial Jordan canonical form. The resaons 
underlying this result are then discussed in view of the results of the paper. 

1.1 CFT 

The CFT corresponding to critical dense polymers has central charge c = —2 and is a 
logarithmic CFT. It is the first member £(1,2) of the logarithmic minimal models £{p,p') 
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[4] with central charges 



c=i.^(p-p'y (1.1) 

With respect to the Virasoro conformal symmetry, it admits an infinite number of repre- 
sentations. In general, these representations are not irreducible — some are reducible yet 
indecomposable (we follow here [6] for the classification of Kac representations). 
From the lattice, a representation Vr,s (also denoted shortly as (r, s)), which we call a Kac 
representation, arises for every pair of integer Kac labels r, s in the first quadrant of the 
infinitely extended Kac table 1, whose conformal weights are given by: 

(2r - s)2 - 1 ^ ^ 

A... = ^ 1.2 
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Figure 1. Kac table of critical dense polymers. 
The conformal character of the Kac representation (r, s) is given by 



ri{q) 



r]{q) 



;i.3) 



corresponding to the Virasoro character of the quotient module K,s/K,-s of the two highest- 
weight Verma modules Vr^s = VAr,s = ^Ar.s- 

These characters are obtained in the limit as ^ oo from finitized characters [5] 



-c/24:+Ar,s 



N 

{N~s+r)/2 



N 

{N~s~r)/2 



;i-4) 



where [^]^ is a g-binomial or Gaussian polynomial. 

A priori, a Kac representation can be either irreducible or reducible. 

Among these are the irreducible Kac representations 

{ir,kp'),{kp,s); r = l,2,...,p; s = 1, 2, . . . k E N} (1.5) 

Since their charactersall correspond to irreducible Virasoro characters, these Kac represen- 
tations must indeed themselves be irreducible. The set (1.5) constitutes an exhaustive list 
of irreducible Kac representations. 

Two Kac representations are naturally identified if they have identical conformal weights 
and are both irreducible. The relations 

{kp,p') = {p,kp') (1.6) 

are the only such identifications. 

We also encounter fully reducible Kac representations 

{{kp.k'p'); k,k'eN + l} (1.7) 

Finally, the Kac representations 

{(ro,So); ro = l,2,...,p-l; = 1, 2, . . . ,p' - 1} (1.8) 
are reducible yet indecomposable representations of rank 1. 

2 Critical Dense Polymers 

We will consider in this paper an exactly solvable model of critical dense polymers on a 
square lattice [5]. The degrees of freedom are localized on elementary faces, which can be 
found in one of the following two configurations: 



(2.1) 

where the arcs represent segments of the polymer. The elementary faces belong to the planar 
Temperley-Lieb algebra [3], and therefore satisfy the following simple equations: 




where the dashed lines indicate that the corners and associated incident edges are identified. 
The parameter (3 represents the loop fugacity which, for critical dense polymers, is set to 
zero. This means that the polymer is not allowed to form closed loops. Therefore it passes 
twice through each face of the lattice, and in the continuum scaling limit it is dense or space 
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filling, in the sense that its fractal dimension is 2. 

The transfer matrix is built out of local face operators or 2-tangles X{u) and boundary 
1 -triangles. 

The local face operators are defined diagrammatically in the planar TL algebra: 



X(u) 



cos(n) 



sin(n) 



(2.3) 



which means that the weights assigned to the elementary face configurations are 



W\ 



cos M 



sm M 



(2.4) 



The local face operators satisfy the Yang-Baxter equation as well as an Inversion Identity. 
The (1, s) boundary 1-triangles are defined as the following solutions to the boundary Yang 
Baxter equation [4]: 




s — 1 columns 



The YBEs, supplemented by additional local relations, are sufficient to imply commuting 
transfer matrices and integrability. To work on a strip with fixed boundary conditions on the 
right and left, we need to work with column Double- row Transfer Matrices represented 
schematically in the planar TL algebra by the iV-tangle 




Where A = | is called the crossing parameter. 

This schematic representation in the planar TL algebra needs to be interpreted appropriately 
to write D{u) in terms of the generators of the linear TL algebra and to write down its 
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associated matrix: 



D(u) 




J= -1 1 



or algebraically: 



N-l 



N-1 N 



j=0 



j=N-l 



being 



Xj{u) = Ij cos(u) + sin(u) 



(2.7) 



(2i 



(2.9) 



(2.10) 



e. = ^ (2.11) 

where e„i is an auxiliary generator which proves useful to express the transfer matrix in 
the linear Temperley Lieb Algebra, the loop generated by such auxiliary generator when 
squaring the transfer matrix, has vanishing a value and one has to remove it by hand in this 
representation. 

Clearly the Sj satisfy in the linear Temepley Lieb Algebra : 



(2.12) 
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(2.13) 



Moreover the boundary tangle Kr^s{u-,C) 2.5 is proportional to the identity in the linear 
TL algebra for (l,s) boundary conditions and does not depend neither from u nor from 
the column inhomogeneities the reader interested to generic (r, s) boundary conditions is 
referred to [4, 29]. 

The matrix representation of the N-tangle is obtained by acting from below (or above) on a 
basis of link states with s — 1 defects, for example the following represents a link state with 
11 nodes and three defects 



(2.14) 



(r',s') = (l,l) 



(r,s) = (l,3) 
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Figure 2. A typical configuration on the strip showing connectivities. The action on the link 
state is explained in the next section. The boundary condition is of type (r', s') = (1, 1) on 
the left and type (r, s) = (1, 3) on the right so there are £ = s — 1 = 2 defects in the bulk. 
The strings propagating along the right boundary are spectators connected to the defects. 



For (1, s) boundary conditions the transfer matrix acts on link states with i = s — 1 
defects which have to be glued into the (1, s) boundary triangle as exeplified in figure 2. 

2.1 Inversion Identities 

For (l,s) boundary conditions the transfer matrix satisfies an inversion identity [5], which 
by virtue of commutativity is satisfied also by its eigenvalues: 



D{u)D{u+-) = [ j-^ r^-Y) =nu) (2.15) 

2 \ cos^ M — sm (u) 
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Such an identity does not depend on s and it can be solved exactly for finite A^, yielding a 
number of solutions which is larger than the size of the D matrix. 

The idea behind the solution is the observation that J-'{u) is an entire function of u which 
can be factorized exactly. The eigenvalues D{u) are determined by sharing out the zeros of 

between the two factors on the righthand side of (2.15). 
The function JF, due to being a square, has only double zeroes which we can define through: 

J'ij + tVk)=0 (2.16) 

where 

i;fc = -^logtan(|) (2.17) 
being tj = ^ for even N whereas tj = for odd A^. 

It follows then that the factorized form of the eigenvalues is [5] for even = 2L: 

tt/c ttJv 

D{u) = 2L2^-2^||(cosec(— ) + efcsin(2M))(cosec(— ) +;Ufcsin(2n)) (2.18) 

k=l 

whereas for odd = 2L + 1: 

TT- vr 2A; — 1 n 2k — 1 

D{u) = 2-2^]^(cosec(---— ) + efcsin(2n))(cosec(-— — ) + /i,sin(2^z)) (2.19) 

k=i 

where 



^li = ei = l (2.20) 

such solutions, however, are too many and one needs to impose some selection rules to pick 
the correct (1, s) conformal boundary conditions. 

The different sectors are chosen by applying selection rules to the combinatorics of zeros. 
A more detailed description of how to choose the e^^^ik in order to obtain (l,s) boundary 
conditions is given in the following section. 
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A typical pattern of zeros for the eigenvalues for = 12 is: 
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A single zero is indicated by a grey dot while a double zero is indicated by a black dot. 



2.2 Selection Rules 

A two column configuration is a couple {l\r) of vectors of length m,n respectively with 
integral entries arranged in decreasing order. 

A two column configuration is called admissible if, calling m the length of I one has: 

lk<rk, k = l,...,m (2.22) 

Notice that configurations with n < m give rise to diagrams which are not allowed. 
It follows then that to each zero pattern of the eigenvalues it is possible to associate only 
one admissible two-column configuration that can be described as in figure 2.23, where one 
is describing the state (3|4,3, 1) . This sequence of numbers is obtained by assigning an 
integer height to the occupied sites in the left column of the diagram, and similarly the rj 
are used as label for the occupied sites in the right column. 
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(2.23) 

The label k in lk,rk is the same k as in 2.18 (we will understand this better from the lOM) 
and one has: 

e„ = — 1, a n E {/i, . . . , Im}, = 1 otherwise (2.24) 

= —1, if n G {ri, . . . ,r„}, /i^ = 1 otherwise (2.25) 

We recall from [5] that the set is the set of all admissible two column diagrams of 

height M with m occupied sites on the left and n occupied sites on the right. 
To each two column diagram V is associated a weight: 

;(P) = 5^/, + 5^r,- (2.26) 



one then defines: 
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one then has the following Fermionic formuals for the finitized characters [5]. 
For odd s one has: 
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For even s, one has: 
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[2.27) 
:2.28) 



[2.29) 



(2.30) 



is the character associated to the set A^^ with respect to the weight 



introduced above. 

From these expressions one can read off at first sight which two column diagrams are allowed 
to contribute to a given sector. The above characters can be reduced to the form 1.4 by 
means of the identity 



M 
m , n 



im(m+l)+in.(n+l) 
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1 L 
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M 
n+1 



M 
m—l 



(2.31) 



Again, for more details con characters and selection rules we refer the reader to the original 
work [5]. 
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3 TBA and Integrals of Motion 
3.1 Derivation of TBA 

The functional equation for the eigenvalues of critical dense polymers is 

^/ TT, /cos^^fu) -sin2^(M)\2 ^ ^ ^ 

D{u)D{u+-) = ( ^ rr-Y) =nn) (3.1) 

2, V cos^(uj — sm [u) / 

the derivation of TBA equations follows closely the work on Ising [23]. The difference being 
essentially that since we are at criticality one has to use fourier integrals instead of fourier 
series. 

First of all we define: 

u = \+'-x (3.2) 



Di{x):=D{u) (3.3) 
^i(x):=^(z|) (3.4) 
one then has that the inversion identity takes the form: 

Di{x - i^)Di{x + i^) = J^i{x) (3.5) 

The function -Di(x) has real zeros and we shall use auxiliary functions to remove the un- 
wanted zeroes: 

p(x, Vk) = i tan(-(x - 2^;^)) (3.6) 

which satisfy 

p{x + i-,Vk)p{x - i-,Vk) = I (3.7) 
one then factors the zeroes in the following way: 

Di{x) := Danz{x) JJ p{x, Vk)p{x, -Vk) (3.8) 
fcex> 

where V = (/i, . . . , /m|ri, . . . , r„) is an allowed two column diagram, whereas Danz is an- 
alytic and non zero (ANZ) in the domain |Im(x)| < |. Danz then by virtue of 3.7 still 
satisfies the same functional equation: 

Danz{x - i-)DANz{x + i-) = J^i{x) (3.9) 

One then Fourier transforms the logarithmic derivative of 3.9 and proceeds precisely as in 
[27, 23] to obtain the TBA equations: 

logL'i(x) = '^\og{p{x,Vk)p{x, -Vk)) + k *\ogTi (3.10) 
kev 
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with k being the usual convolution kernel 



k(x) = ^— - (3.11) 

^ ^ 27rcosh(a;) ^ ^ 

such a kernel arises from an integral of the type: 

» iakx 1 

dk^^ w = -7^ 1 T (3-12) 

e/3fc + e-/3fc 2/5 cosh (^) ^ ' 

which is evaluated using the residues method, and actually can be used as a tool to figure 
out how the fourier tranform is defined. 

We now recall the definition of continuum limit for this model, this can be obtained by 
simultaneously going into the braid limit u —>■ ioo and in the thermodynamic limit N ^ oo 
and simultaneously keeping the following quantity fixed: 

Ne'^''' = const (3.13) 

this is obtained if u diverges as u ~ | log A^. 

To obtain a meaningful expression in the continuum limit one has to subtract the divergent 
part out of jFi. so one has for even N 

T{u + - logiV) (cos(4e2») - 1) (3.14) 

whereas for odd 

^(n + -logiV) (cos(4e2») + 1) (3.15) 

the scaling limit of the p functions is denoted by p and is computed by using the exact result 
for the 1-strings: 

ffc = "2 logtan(^) (3.16) 

being tj = ^ for even N whereas tj = ^^■^^"^''^ for odd A^. It turns out that for even one 
has: 

1 AjTt 

p(x, -Vk) = tanh(-(x + log (^))) (3.17) 

and 

p{x,Vk) = l (3.18) 
the convolution term after subtracting the explicit A^— divergent term looks like: 

k*Ti^ / log(cosh(4e-^) ± 1) (3.19) 

i_oo 27rcosh(x-t/) 

this convolution product however still hides a divergence which we want to remove. 
For this reason we regulate the convolution term with a cut off: 

r^'^ dv 

k*T^= / -f log(cosh(4e-^) ± 1) (3.20) 

J_iog(^) 27rcosh(x - y) 

12 



we then expand the hyperbohc cosine, change variables and integrate twice by parts: 

* ~ y (_l)n+lg(2n-l)x / ^g-(2n.l), log(cosh(4e-^) ± 1) = 

= E(-l)"^';^^^'""'^^ / ^^'""'^ log(cosh(t) ± 1) = 



n=l 

°o 1 / 1 /-Af j-2n 

^} ' 7r4(2«) \ 2n{2n -I) Jo cosh(t) ± 1 



{3.21] 



n=l 



^ .N^n ""'^^J^ + ^JL^Ar2n-i iog(cosh(iV) ± 1) 



2n(2n - 1) cosh iV ± 1 2n - 1 
we now define: 

logZ±(x,iV) = e(2«-i)-( -i -iV^" 'i''^^ +^^iv2»^ilog(cosh(iV)±l)) 

^ ^ ' ^ ^ 7r4(2") ^ 2n 2n-l coshA^il 2n - 1 v ; 

n=l ^ ' 

(3.22) 

which can be summed up as: 
logZ+(x,A^) = 

_^21og(l + j^e^^N'^) sinh(iV) + arctan(e^'f )((1 + cosh(iV)) log(l + cosh(iV)) - iVsinh(iV)) 
^ 47r(l + cosh(A^)) 

(3.23) 

and 

\ogZ~{x,N) = -^f2e-^coth(^)log(l+-^e2^Ar2)+arctan(e^^)(-iVcoth(^)+log(cosh(iV)-l)) 
47r V 2 16 4 2 

(3.24) 

these functions Z^{x, N) provide us the divergent part we want to subtract, for this reason 
we define: 

^finite 



\ogDf,r^ue{x,N) = \ogD^{x,N) - \ogZ^{x,N) (3.25) 



so that by calhng D{x) the continuum scaled version of DfinUeix) one has the following 
continuum TBA: 



dt- 



„=i ■■- 2n{2n-l)Jo cosh(t) ± 1 

(3.26) 

in this form the convolution kernel is hidden, it could be recovered by further resumming 
the series, but in this form it turns out that it is ready for use in the next section. 
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3.2 Integrals of motion 

We shall now deal with the expansion of D which will yield the eigenvalues of the BLZ 
involutive charges. Again, as in [23] [25] we consider the following expansion: 

oo 

log^(x) = - 5^ f/„/2„,„ie(2-i)- (3.27) 

n=l 

where the /2n-i are the eigenvalues of the BLZ involutive charges. 
For the auxiliary functions one uses an expansion like: 



\ (2n-l) 

logtanh(^i^) =m-2^e(2"-i)-ili- (3.28) 

n=l 

We must be careful, and consider both the contributions of single and double 1-strings, in 
the case of double 1-strings the summation term carries an additional 2 factor coming from 
the log of a square. 

The convolution term, which has already been expanded in the previous section, can be 
further simplified by means of the following identies (see appendix): 

^^TT TTT = 4^(1 - 2'"'")r(2n)C(2n) (3.29) 

1 + cosh(sj 

and also, in the other case that the useful integral is: 

oo 2n 

ds — — = 4nr(2n)C(2n) (3.30) 

cosh(s) - 1 V ^■'v / 

So that one now reads off in one case: 

UJZ-i = ^42n-i2„(2n - 1) X "^"l + cosM.) ^^'^^^ 
whereas on the other hand: 

UnP2n-l = ^42n-l2n{2n - 1) J, '^'cOsh(s)-l ^^'^^^ 

in view of the lattice selection rules explained in the previous section , it will prove useful 
to rearrange the 1-strings according to the double column diagram description. From this 
point of view there is actually no difference between lattice selection rules and continuum 
ones, and actually the j will be nothing but the labels of the two column diagrams. 
So that piecing up highest weight and excitations one gets for the —1/8, 3/8, .. . family: 



C/„/,„.. ^ 2(^) '^"-"(5^1: (2,-1) '^"-\(-l)"(l-2'-.)r(2.„-l)^) (3^33) 
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whereas for the 0, 1, 3, . . . family 

UJ.„-. ^ 2(9 ''"-'(^^ (2,)-" . i-imn - D^M) (3.34) 

j&V 



Notice that an immediate check can be made by considering the case n = 1 which corresponds 
to the energy Ii of conformal field theory, indeed it is not too difficult to realize that f/i = vr 
gives for the energy of a state described by the diagram V G Am,m+2k in the odd A^— parity 
sector: 

h{V) = (w{V) - m - A; - ^ + 1) (3.35) 
whereas the same expression for even A^— parity is: 

I,(V) = (wiV) + ^) (3.36) 

these expressions explicitly exhibit the correct value of the central charge and the related 

conformal weights 0, — |, and thus are a good hint of what is going on. 

So by applying the (1, s) selection rules to the quantum numbers one obtains all the correct 

CFT characters, we shall see by the way that continuum selection rules are obtained from 

the lattice selection rules by allowing for infinite height two-column diagrams. One then has 

that comparing with explicit diagonalization of the matrices obtained from appendix 9.2 for 

the BLZ lOM, that the above formulas describe for generic n the eigenvalues of the local 

involutive charges of BLZ, provided that we isolate the value of the constants f/„. 

The above expression is identical to Ising, aside from j taking values into fermionc partitions 

V which are essentially one-column diagrams. 

The first few constants are found by direct comparison with CFT to be: 

f/i = TT (3.37) 

U2 = ^ (3.38) 

f/3 - ^ (3.39) 

Actually one notices that these constants are precisely those one can obtain from the Ising 
model by requiring to describe the c = —2 theory instead (aside from the factor of 2 appear- 
ing in front of the formula). Actually the above expression is the same expression as in the 
Ising case, this is true of course if we consider the largest eigenvalue. The excitations are a 
bit different due to the presence of double zeroes. 

It may be instructive to observe that the behaviour of excitations in this model is actually 
encoded in some properties of the Bernoulli polynomials. If we take the BLZ formula for 
highest weight eigenvalues in the first column of Kac's table: 

IZ^, = 2""i?2„(^) (3.40) 
and use the following property of Bernoulli polynomials 

B2n{-^) = i?2n(^) + (3-41) 
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we realize that i?2n(0), -B2n(|) (so s = 1,2) represent the highest weights 0, — |, and that 
going beyond those in Kac's table implies adding a certain number of 1-strings which is at 
this point trivially guessed. 
Suppose for example s to be odd, one then has: 

= 2'-"n ^ + 2-"52„ (3.42) 
i=i 

Actually this shows also that the odd power behaviour of 1-string contributions is actually 
encoded in the BLZ formula and therefore this formula alone should be enough to suggest 
the structure of all the excited states. 

One can be even more explicit, and resum the contribution of the quantum numbers, to get 
the explicit expression for Ai^^ and its close relatives pertaining to the higher lOM: 

Tvac ^ (s-l)^(s-3)^ l_ 

3 64 120 
Tvac ^ {s-lf{s-2f{s~3f , J_ 
5 512 ~^ 336 
Tvac _ (s-l)^(s~3)^(ll-8s+50s^-24g3_)_334) ^ 

7 ~ 12288 480 



and, in general: 

IZ 1 = 2^-"n y (-l)'f2n^^ (2^-1)'- (, _ 3)fc (3.44) 



k=0 

The constants Un do not depend on s and can be fixed either from the highest weight 
Ai 1 = or from the coefficients of 1-strings in the expressions obtained from TBA: 

This simple observation allows to obtain immediately the closed form for the coefficients Un- 

f/„ = — — (3.46) 

2"-i(2n2-n) ^ ^ 

We shall see that the constants f/„ will reappear in many different places in the following 
sections, and thus their explicit knowledge is an extremely important basis for all the sub- 
sequent analisys. 



3.3 Euler-Maclaurin and Integrals of Motion 

The goal of this section is to extend the Euler Maclaurin analisys carried out for polymers 
in [5] and for Ising in [21] to all orders on 1/A^. As we already know the first order turns 
out by general arguments to be proportional to the eigenvalues of the energy Ji of the 
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underlying CFT, and is used as a tool to identify the central charge and to prove that the 
finitized characters yield for N ^ oo the quasi rational characters of CFT, after extracting 
the divergent and constant parts which are proportional to the following bulk and boundary 
free energies: 

fbuikiu) = log V2 - - / ' log(-^ + sin 2u)dt (3.47) 
TT Jo sm t 

fbdy = log{l + sm2u) (3.48) 

The higher order corrections turn out to be related to the conserved quantities of CFT as 
well by the following asymptotic expansion: 

oo ^ 

log D{u) = -2Nhuik{u) - fbdyiu) + sin(2M)P„(sin(2M))f/„J2„-i (3.49) 

n=l 

where P„ are polynomials whose explicit form is: 

n 

Pn{a) = 5^(-l)("+^)C„,fca2('=-^) (3.50) 

k=l 

being 

and the succession 6„ is: 

2(2n-l) 

This result is remarkably compact and independent of the parity of A^, and actually up to 
first order it is well known to be a general feature of RSOS models and Logarithmic Minimal 
Models. The somewhat surprising simple result obtained here to all orders relies a lot on the 
factorization of the eigenvalues, and it would take a lot of additional effort just to investigate 
the persistence of such a property in the general case. 

We will now give a brief explanation of how this calculation proceeds in the case of even N . 
First of all introduce the auxiliary function F, defined as: 

Fit) = log(tcosec(t) + t sin(2M)) (3.53) 

in terms of this function one can express the logarithm of the eigenvalues in the following 
form: 



L-l , L-l , 

\ogD{u) = (1 - 2L) log 2 + log(2L) + 2 F(— )) - 2 log(— )) + 

k=l k=l 

v-j 1 - sin(2n) sin(|f ) 1 - sin(2n) sin(|f) 

h 1 + ^i^(2n) sin(t) 1 + sin(2n) sin(||) ^ 



(3.54) 
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the sum over F is evaluated by means of the Euler Maclaurin formula: 



(3.55) 

the sum over the logarithms is evaluated by using the asymptotic of the F function: 

E log(|^) 4 log L . L(log f - 1) . log 2 H- f (^)-.^-J?^(^)-. (3.56) 

fc=l n=l ^ ^ 

one then uses the values of the derivatives of F 

F(|) = log(|) + log(l + sin(2M)) (3.57) 



.2. iwTT, 2(2fc-i)(2(A;- 1))! 



^(2fc-i)^Q^ = sin(2M)Pfc(sin(2M)) (3.58) 
,7r _ 2(^fe-^)(2(A;- 

and notices that F*^^^~^)(^) is engineered to cancel the contribution of the T function, whereas 
the even derivatives of F drop out of the calculations regardless of their explicit form. 
The excitations are included by noticing that they are generated by: 

(-0) 

n=l ^ ' 

And finally piecing up one arrives at 3.55. 
If now one calls 

oo ^ 

Gn{u) = J2 sm{2u)P^{sm{2u))Uj2n-i (3.61) 

n=l 

it is possible to reshuffle the sum so as to collect a given power of sin (2m) as: 



Gn{u) = J2MN)sm^^-\2u) (3.62) 



1=1 



^■'W = f ""'-tz-'r (3.63) 



being 

oo ^ 

JY{2r-l) 
r=l 

actually one can do even more, and resum the above series explicitly. 

The expressions one obtains essentially depend on the parity of A^, which is conveniently 
parametrized for even N as N = 2D + 2, whereas for odd N as N = 2D + 1. 
It is also convenient to isolate the constant and divergent contributions, as well as the 
contribution of the excited states: 

Ki{D) = Ki{D) - Kf'' (3.64) 
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Ki{D) = Kr\D) + Kf\D) (3.65) 
the K'^^^ and K^^'^ are defined independently of tlie parity of N: 

Tin-\)V{n) 1 

™= i2,.-\)2^^-^ ^ i:(-ir^'(!:;^)-((2™+i)<,) (3.67) 

whereas the other pieces are, for even TV: 



Kin) ^ ^.-^l^ gs. (i^) eosee((2..1)|) s.n 

(3.69) 



m=0 

whereas for odd TV one has: 



{2j - l)vr {2j - l)n 

2N ~2{2D + 1) ^^-^^^ 

KiD) ^ ^ + (2.-!)2^n-j /-^r" C:;.') -ec((2^+ l)tO (3.71) 

In the next section we will recognize the as eigenvalues of suitable iV— tangles defined 
in the Temperley Lieb algebra. 

We now want to resum the contribution of the divergent part, for reasons that will become 
clear in a short time: 



sin--(2n) = 1 log {^-^!^) + ^-^iM2u) + l ((1, | 2); (| \), ^n\2u)) sin^(2n)) 

n=l + I J ^ 

(3.72) 

where 



is the generalized hypergeometric function. 

It is indeed remarkable that the bulk and boundary free energy produce very neat cancella- 
tions with the resummed divergent part, by means of the following identity: 



dt logfcosec it) + sin(2n)) = - log(l + I - sin^(2n))- 

2 ^ (3.74) 

+ sin(2n) + ^ 3i^2 ((1, | 2); (| ^); sin2(2n)) sin3(2n) 
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one then uses also the following expansion 

OO I 

log(l + ^1^) = ^^, z'- (3.75) 

and ends up with the following expression for the eigenvalues D: 

^ A 

log = 5^ ^ sin" (2m) (3.76) 



where 



^.„^(2„).|i-H-iV(-ir'^-^^p^|^| (3.77) 



/l2„_i = (2n - l)!if„ (3.78) 
one then introduces the complete Bell polynomials: 

n=0 

which are defined recursively as: 

n 

B„+i(Ai,...,A„) = (^) A„-fc+iSfc(Ai,...,Afc) , Bo = l (3.80) 

fe=0 

^3 + 3A1A2 + A3 



54 =A1 + ^A{A2 + 3A^ + 4A1A3 + 

^5 + lOA^A^ + ISAiA^ + 10^2^3 + 5Aiv44 + As 



n=0 



(3.81) 



in terms of these polynomials one has the following expansion for the eigenvalues: 

D{u) = ^"^^''V"'^"^ '^^"^2m) (3.82) 



actually it is possible to read off from the factorized form of the eigenvalues that they are 
polynomials in the variable sin(2u), whereas the above expansion is an infinite series. This 
is due to Euler Maclaurin (which was our starting point) being an asymptotic expansion. 
Fortunately this is not a problem. It turns out that one simply has to truncate the above 
expansion to get the exact result: 



2D 



n=l 
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This decomposition will be lifted from the eigenvalues to the transfer matrix itself in the 
next section. 

It is also worth, again for the meaning it will carry in the next section, to recast the inversion 
identity in the following form: 

BJ2^A2,...,2j§^A2k) 
D{u)D{u + X) = Y, ^ — — -sm^\2u) (3.84) 

while we are about it we also give the following explicit evaluation of the above Bell poly- 
nomials, which can be obtained by explicitly expanding jF(n): 

B, (2-A2, . . . , 2^A2fc) = k\F2k (3.85) 



being 



2D 



F2m — fr,mgr,N (3.86) 

r=0 

f V (-l)-^^r(l + l) .rx 

A,m-2.2.r(m + i)r(i + i-m) U) ^^•'''^ 

r 

3r,N = hm,Nhr-m,N (3.88) 

m=0 

Acyl I N 



''m,N 



((-1)^ - 1)2'"-^ m>N 



(3.89) 



4 Integrals of Motion on the Lattice 

In this section we want to put the attention on the meaning of those misterious results which 
we obtained from Euler Maclaurin. 

What happens is that the transfer matrix admits the following expansion: 

D(n) = l + y BrA^^^^^y^^n) g.^„^2n) (4.1) 

n=l 

where, following the notation of the previous section we define 

A2„ = A2nl (4.2) 

A2n-i = (2r2-l)!K„ (4.3) 

the and K„ rightfully deserve to be called Lattice Integrals of Motion, and they are in 
involution: 

[Kz,K,„] = (4.4) 
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by construction they are diagonal in the same basis as the transfer matrix itself, so that if 
we label an eigenstate by the corresponding 2-column diagram X> one has: 

K^\V) = Kn{V)\V) (4.5) 

B{u)\V) = D{u)\V) (4.6) 

where Kn is the quantity which we computed in the previous section. 
A comment is in order to clarify why the K„ deserve to be called lattice integrals of motion, 
this is because their eigenvalues which we computed from euler maclaurin analysis in the 
previous section behave to leading order as the local involutive charges of the underlying 
conformal field theory: 

Kn{V) ~ ^j^hn-i{V) + 0{N-'-') (4.7) 

where all the constants of proportionality appearing in the above aymptotics have exact 
expressions reported in the previous sections, and hence the importance of all that otherwise 
meaningless calculations done with Euler Maclaurin. 

We notice also that at this stage the constants f/„ 3.46 obtained from the continuum TBA 
and expansion into BLZ local integrals of motion reappear in this simple formula and hence 
once more we recall the importance of their explicit knowledge. 
The tangles K„ and K„ are related by subtraction of a diagonal divergent part: 

Kz = - Kfn (4.8) 

We are now going to exhibit explicitly how the lattice lOM are built from the generators of 
the TL algebra. 

First of all we introduce the boundary symmetric TV— tangles B^: 

Bfc = Cfc + eN_k (4.9) 
and the following nested commutators, which for even TV = 2D + 2 take the form: 

2D+3-2n 

Hn = ^ [ej, [e^+i, [ej+2, [• • • , [ej+2n-3, ej+2n-2] • • •]]]] (4.10) 

while for odd = 2D + 1 the bound of the summation is different: 

2i:>+2-2r!. 

H„ = ^ [ej, [e^+i, [ej+2, [■ ■ ■ , [ej+2n-3, ej+2n-2] • • •]]]] (4.11) 
i=i 

The idea of introducing nested commutators in TL expansions is not completely new, for 
example it has been used in [28]. 

In terms of the H„ one has the following form for the fist few lOM: 

Ki = Hi (4.12) 
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K2 = — H2 + -Ki - — Bi (4.13) 

12 6 12 ^ ' 

K3 = — H3 + — H2 - — [Ki, [Ki, Bill + — Ki - — Bi - — B2 (4.14) 

^ 80 20 80^ ' 40 80 80 ^ ^ 

K4 = H4 + H3 + H9 fK2, [Ki,Bill [Ki, fKi,Bill [Ki, [Ki,B2llH 

448 224 448 ^ 112^ ^ 64^ ^ " 224^ ^ 

:Ki — TT^Bi — — B2 — TT^Ba 



112 448 28 448 

(4.15) 

K5 = H5 + H4 + Hs + H2 [K2, [Ki, Bill fK2, [Ki, B2II + 

^ 2304 ^ 288 576 ^ 288 192^ ' 96^ ' 

- ^[Ki, [Ki,Bill - ^[Ki, [Ki,B2ll - ^[Ki, [Ki.Bgll - ^[Kg, [Ki,Bill + 

35 — 1 11 19 ^ 7 ^ 

H Ki H Bi Bo — Bt — B4 

1152 128 384 1152 2304 

(4.16) 

We consider also the inverse relations which give the nested commutators in terms of the 
boundary tangles and lOM: 

H2 = I2K2 - 2Ki + Bi (4.17) 



H3 = 8OK3 - 48K2 + 2Ki + [Ki, [Ki, Bill - 2Bi + 3B2 (4.18) 
H4 = 448K4-480K3+108K2-2Ki+[12K3+Ki+Bi, [Ki, Bill + [2Ki, [Ki, B2II-4B1-2B2+5B3 

_____ ___ _(4.19) 

H5 =2304K5 - 3584K4 + I6OOK3 - I92K2 + 2Ki + [-8B1 + 5Ki - 36K2 + 8OK3, [Ki, B1II + 

+ [-Ki + 24K2, [Ki, B2II + [3Ki, [Ki, B3II + 14Bi - 2B2 - 2B3 + 7B4 

(4.20) 

this structure shows some remarkable regularities, indeed it is possible to suggest that the 
general structure should be something of the form: 

n n—1 n—2 

Hn = J2^l,nKl + J2^l,nBl + J^lPlA^^^ ■ ■ ■ ^^n-l-l, Bi), [Ki, B^H (4.21) 
1=1 1=1 1=1 

where 

n-l~l 

Pi,n(Ki, . . . ,Kn-i-i; Bi) = J2 Pi,n,hKh + ai^nBi (4.22) 



h=l 



in particular 



and one has also: 



P„_2,„(Ki) = {n- 2)Ki (4.23) 
Sn-l,n = 2/2-3 (4.24) 



c,„=(-i)'-(2/-i)2--i i (4.25) 
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so that the task of solving the problem is reduced to identifying the Si^n, Pi,n,h, CLi,n successions 
appearing in 4.21. 

It is also worth considering the inversion identity in the Bell polynomial form: 



^bJ2'4A„...,2^A,,) 
1 + > ^ r-. — si 



D(«)D(m + A) = 1 + ^ ^- — — ^-sm^\2u) (4.26) 



A;! 

k=l 



one notices that on the right hand side of the above equation only the even can appear. 
What turns this identity into an inversion identity is simply the fact that the are pro- 
portional to the identity. 



4.1 A Different Point of View 

In the previous section we have dealt with representations of the lattice lOM in terms of 
nested commutators, and we got to the point of giving an ansatz for their general form. This 
description, however does not capture the full picture, since for some fixed size D there is 
only a finite number of independent lattice lOM, the higher ones being dependent from the 
lower ones from some point on. We shall understand better the degeneration of the lOM 
from a different point of view. 

Therefore we give a name to something we have already introduced: 

- i?„(Ai , Ao, . . . , A„) 

D„ = ^' , '—^ 4.27 

n! 

where the identity: 

2D 

D(m) = 1 + ^ D„ sin"(2M) (4.28) 

n=l 

can be explicitly proven by first starting from the more natural expansion, which is readily 
obtained from the expansion of the elementary faces in terms of connecions: 

^ 2N-1 

D(n) = -J2 cos2^-'=-i(u) sin'^-^(n)Dfc (4.29) 

k=l 

Now, use of crossing symmetry and some relabelling of summations tells us that the D„ are 
related to the D„ by^: 

I - 1 

6. - . ± '-^'''"7-'"^^' (" D.«-.} (4.30, 

^the author thanks Jorgen Rasmussen for proving a part of the foUowing identities 
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and 

N-l 

At=l d,a,P>0 

ei...e^+d-i[ ^i---^ic){ ^Jp---^n)}+ (4.31) 

li+d<ii<...<ia<N-l IJ'<jl<---<ji3<N~l 
li<ii<...<ia<N~l ^i+d<ji<...<ji3<N~l 

The expressions for the first D„ read: 
Di = Bo 

D3 = Bo + 4 Y {e^'^i} 

l<i<j<N-l 

D4 = 4Bi + 8 5^e, + 4 Y {e^, {e,-, 6^}} (4-32) 

j=2 l<i<j<A:<Ar-l 

D5 = Bo + 4 Y {ei,Bi} + 8 ^ {e„ej} + 4 Y 

2<i<N-2 '2<i<j<N-2 l<i<2j<N-l 

+ 4 Y {e»,{ej,{efe,ei}}} 

l<i<i<fc<«<Af-l 



where 

a <n b <^=^ b — a > n (4.33) 
Let us introduce the inverse Bell Polynomials, which are defined by: 

A:=l fc=l 



defined by recurrence as: 



^n+l(C*l, . . . 5 C*n+l) — Cn+l ~ -j^ ^ C*n-A;+l^A;(C'l, . . . , C**;) , ^0 — 1 (4.35) 

by calling 

Cfc = klBk (4.36) 



one has by definition that 

_ Y2n-1 

~ (2n- 1)! 
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and 



2n ^ ^ 2n2r(i-n)r(-n)r(r2+l)y (2n)! 

the second relation allows one to eliminate D2fc from the definition of the involutive charges 
by solving for the linear term. 
One then has: 

Ki = Di 

3!K2 = 6D3-D? + ^(iV-2)Di 

5!K3 = I2OD5 - eOD^Dg + 30(iV - 2)D3 + 9Df - 15(iV - 2)D? + 15(— - j - l)Di 

7!K4 = 504OD7 - 2520D^D5 + 1260(iV - 2)65 - 2520DiD^ + ISQOD^Ds - 1890(iV - 2)D^D3+ 

+ 315(^ - y - 2)D3 - 225DI + ^(iV - 2)D? _ ^(iV^ _ 3iV + |)D?+ 
105 

+ (iV^ + 3iV2 - lOA^ - 24)Di 



(4.39) 

from these expressions it is now obvious that for fixed D = for k > 2D and the 
expressions for the involutive charges have to degenerate correspondigly. Oviously this fact 
is not trasparent from their expressions involving commutators and boundary A^— tangles, 
this form for the involutive charges is suitable only for the non exponentiated form of the 
transfer matrix, where only the non degenarate charges do appear, since all higher order 
terms cancel due to the degeneracy pattern. 

One can also chose to express the involutive charges in terms of the A^— tangles, although 
one has to be careful in doing so, because first the degeneration of the above expressions, if 
any, has to be obtained by putting the appropriate to zero, and then one can proceed 
use (4.30). 

In the non degenerate case one has: 



Ki = 


D2 

4 












3!K2 = 


8 












5!K3 = 


8 ^ 






1024 2 8 


D2 




7!K4 = 


315_ 
16 ^ 


315_2^ 
D2D6 

128 ^ 




315_ _2 
= - 128°'°' 


945 
^ 2048 


D^D4 


+ 


315_ 


225 

16384 2 


315_, 
128 ^ 


315_ 







(4.40) 

315_2^ 
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we can also invert these relations and obtain the relations between the involutive charges 
and D„: 

Di = 2 1 

D2 = 4Ki 

D3 = 2 1 + 4Ki 

D4 = I6K2 + 

D5 = 2 1 + 32K1K2 + ^K^ - 4Ki 

De = 64K3 + 32K^K2 - 32K2 + - + 4Ki ^^'^^^ 

15 3 

D7 = 2 1 + 8kJ - 4kJ + - 96K1K2 + ^k'Ks + 64kI + I28K1K3 

Dg = 256K4 + 128kJk3 - 256K3 + I28K1K2 + — K1K2+ 

3 



128K;K2 + 64K2 + — k' - -k! + -K' 

^ 315 ^ 15 ^ 3 ^ 



one can notice that in this representation the coefficients are independent of A^. It is 
also worth commenting that although the variables D„ have the appealing feature that in 
the expansion no high powers of Di occur, which would generate diagrams not allowed to 
appear in a double row A^— tangle, such diagrams will eventually appear in the involutive 
charges and disappear after building the Bell polynomials out of them. On the other hand, 
the commutator representation of the involutive charges does not suffer of the presence of 
high powers of the hamiltonian, but instead is made in such a way that the Bell polyno- 
mials will explicitly contain the unwanted terms. Such terms of course will always end up 
disappearing in the final result for the transfer matrix. 



5 Symplectic Fermions 
5.1 Generalities 

We now want to discuss the different ways of describing states for critical dense polymers 
and their relation to symplectic fermion states in the continuum limit. 
The CFT describing symplectic fermions is built from the following stress energy tensor 
[12, 13]: 

nz) = ^:x{z)-x{^)- (5.1) 
where we have introduced the notation 

X-X = C^a,/3X"X^ (5.2) 
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where da^p is the antisymmetric tensor satisfying (i+ ._ = 1. 
the field 



Xiz) = ( ) (5.3) 



X [z 

is a quasi primary field of scaling dimension 1, and by introducing the mode expansion 



x(-) = E l^r (5-4) 



one has that the modes satisfy the following anticommutation relations: 

{x;^.,X^} = rnrt^H-„ (5.5) 
it follows that the Virasoro modes can be expanded in Symplectic Fermion modes: 

■^n — 2 ^ ] • ■ Xn-m • (5-6) 

m 

where the summation is over Z when the modes are considered to be acting on the vacuum 
fl wereas the summation is over Z — | when the action is over the twisted vacuum On 
notices that in the twisted sector there are no fermionic zero modes. 

The energy Lq does not have a diagonal action in the sense that there exists a logarithmic 
partner uj of the vacuum fl such that: 

Lou = n (5.7) 

Lon = (5.8) 

By the way, is we decide to build the module over the vacuum Q, by virtue of 5.8 the 
logarithmic partner does never appear. 

In order to select the Xi,i character it is necessary to require traslational invariance, which 
implies: 

X-i ■Xo^ = (5.9) 

this can be obtained by requiring: 

Xo^ = (5.10) 
this can be interpreted also as a condition on the fermionic states 6, defined as: 

Xo^ = -0 (5.11) 

X'^e^ = d'^'^Q (5.12) 

so that 

Xo^^ = 2xtXoO (5.13) 

implies that 5.10 is equivalent to one of the 2 component of the zero mode annihilating the 
fermionic state 0. 
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The theory has a global s/(2) isospin symmetry, and the free fermion field x transforms as 
a j = I representation of sZ(2)[12]: 



[J , J J — ZJ 






(5.15) 


1 

[Ax^(^)]=±2X^(^) 


(5.16) 




(5.17) 




(5.18) 


J°fi = J^fi = 


(5.19) 


by virtue of this global symmetry the highest weight states will always fall into irreducible 



representations of s/(2) carrying isospin j G |N: 

where the round brackets denote symmetrization over the upper indexes. 

The whole multiplet can be obtained by acting on |j, — with rising operators: 



(5.20) 



k=l k=l 

For example the state with weight Ai^g 
composed by the following four states: 



2j 



xt2X-i^ = \j, 1- 
(5.21) 

6 forms a j = | multiplet of sl{2) which is 



3 3 
2 

3 1 
2'~2 



(+ - -)n 
XUX-2X4^ 



3 1 
2 

3 3 
2'~2 



(5.22) 



, 2) = AxUxzl^ 

) = XZ3XZ2XZi^ 

1), covering all entries in Kac 



in general the states |j, have confomal weight Aj = j{2j 
table with integer conformal weight. 

All the other entries in Kac table can be described by introducing a twisted vacuum /i and 
using fermi modes labelled by half integers: 



I ■ \ (+ 



1 y 

1 I 



2 



(5.23) 



in this case the multiplet has conformal weight Aj = — | + 2j^. 
So that if one picks one of the \j,m^, either twisted or untwisted, it is possible to build a 
W module over it. 

In order to estabilish connections with boundary conditions corresponding to Virasoro rep- 
resentation one introduces the infinite dimensional Clifford algebra ^(2) of [18] generated by 
the commutation relations 5.5, which admits an SI2 action for which the \j,m^ are highest 
weight vectors. This SI2 action commutes with the Virasoro algebra generated by the stress 
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energy tensor and therefore irreducible ^(2) modules decompose as (see again [18] which we 
are almost quoting literally): 

Xi = TT,- ® V2,+i,i (5.24) 

X2= 07r,®V2,+i,2 (5.25) 
where ttj is a 2j ' + 1 dimensional representation of SI2. 

We can select a Kac representation by applying a suitable projection operator to the modules 
Xi,X2 onto the eigenspaces of fixed isospin, this operator acts in the following way 

n2^-.-(Xi) = l,,„®V2,+i,i (5.26) 

n2^-.-(Xi) = l,,„®V2,+i,2 (5.27) 
notice that it makes sense to state that: 

j 

-^j = h,m (5.28) 

m=-j 

Furthermore the modules Xi,X2 admit an SI2 action which is used to classify the Virasoro 
highest weights int multiplets of isospin quantum numbers j, m. 

So if we consider the highest weight vectors of SI2 ® Vir we can build modules over them of 
the form (and similarly Q^^'"^ modules built over the twisted vacuum /i): 

Q2^'"^= {L_„,...L„„Jj,m)} (5.29) 

these modules are related to the projections we have introduced: 

Q2j,m. ^ n2-''"(Xi) (5.30) 

building these modules explicitly in terms of fermi modes one observes that these modules 
are irreducible (being projections of irreducible modules onto one a one dimensional repre- 
sentation of SI2 quantum numbers j,m), this means that we already removed the first null 
vector which appears in the module, and since the embedding structure of modules for this 
value of the central charge is linear, no further new null vectors are allowed to appear. By 
linear embedding of modules we mean that each new singular vector will appear inside the 
module generated by the first singular vector, which however we have already modded out. 
And all we need to use to build a Q^-^'"* module is the commutator between virasoro modes 
and fermi modes: 

[^-n, X J = /X_(n+Z) (5.31) 

To select a sector labelled by (r, s) it will be necessary to identify which combination of 
fermi modes corresponds to the null vector at level rs, for example both (2,1) and (1,5) 
have A = 1, the difference being that by virtue of some selection rules we will be able to 
identify which submodules we have to throw away. These selection rules will be the same 
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as the lattice selection rules. 

Notice finally that by using 5.28 we can express the irreducible W modules as: 

^1=00 2'^'"^ (5-32) 

^2=00 Q'^'"^ (5.33) 

we will comment in the next section on the fact that these modules coincide with the W 
modules of [29], this will be recognizable from character identities, explicitly: 

Vi^i = Xi (5.34) 

Vi^2 = ^2 (5.35) 
We have thus made contact between notations of [18] and [29]. 



5.2 Selection Rules and Characters 

In order to deal with selections rules let us introduce some obvious notation. 

Let be the set of all admissible two column diagrams with m occupied sides on the left 

and n occupied sites on the right, where each diagram has no height restriction. 

One then introduces: 




(5.36) 



and 

oo 

771=0 

The object x^'^'^Kq) is the character of a constant J° subspace of the irreducible SI2 ® Vir 
module which we called Q(2i,7n) jg b^jit on one of the highest weights |j, mj, the choice 

of m is not important here, because all such modules are isomorphic due to sl{2) invariance, 
and can be generated by the action of the operators, for this reason we will suppress the 
label m in Q2i,777 

These characters are simply related to the (1, 2/c + 1) quasi rational characters by virtue of 
the formula: 

2i 

= ^(-l)2^-^xi,2fc+i(g) (5.38) 

k=0 

which can be inverted to yield: 

XiAHi{q) = X^''~'\q) + X^''\q), J e (5.39) 
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notice that 

Ai,4,+i=j(2j-l) (5.40) 

and that one defines also x''~^Kq) = 0. Since the x^^^\q) are well defined characters with 
positive coefficients, 5.39 can be interpreted as meaning that the Virasoro module Vi,4j+i 
admits the following decomposition: 

Vi,4i+i = Q('^"') © (5.41) 

where, again, Q^~^^ = 0. 

Indeed, it is possible to obtain information on generic decompositions for r > 1, by means 
of the following identity: 

k 

Xi+fc,s(g) = ^Xi,s-2k+4p (5.42) 

p=0 

where it is understood: 

Xi,o(g) = (5.43) 
Xi,_.(g) = -Xi,s(g) (5.44) 

which tells us for example that: 

X3,5 = + Xi,5 + Xi,9 = + X^'^ + X^'^ + X^'^ + X^'^ (5.45) 

which implies for example that the following modules admit decompositions such as: 

V2,3 = © Q^'^ © (5.46) 

V2,5 = 2^°^ © Q^'^ © Q^'^ © Q^'^ (5.47) 
V3,5 = Q^°^ © Q^'^ © Q^'^ © Q^'^ © Q^'^ (5.48) 

It follows from simple cancellations of characters that: 

XnM) = X^'"'\<l) (5.50) 

which gives the identification: 

V„,i = (5.51) 

one the focuses on (n, 3): 

Xn,3 = x^"~'^(g) + x^"~'^(g) + x^"Hg) (5.52) 

so that 

V„,3 = © © (5.53) 

the case of (n, 5) gives: 

xn,5 = x^^-'n?) + x^^-'H?) + x^""'^?) + x^"Hg) + X^""''^?) (5.54) 
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corresponding to: 

V„,5 = Q^""'^ © © Q^"-'^ © Q^") © (5.55) 

Where of course all the Q*-"-* with negative n are empty. 

We are now going to deal with the twisted case. 

Following the selection rules defined on the lattice we define: 

oo 

X^'^^ = q~^^~^Y.1~"'UZ2^ (5-56) 



n 1 

Notice that in order not to introduce further notation we are using for 5.56 the same name 
as 5.37, this should not rise any confusion, since we are working in a different sector of the 
theory. 

One then has that, as in the previous case x^^^ = for n < 0. On the other hand, whenever 
the x*^""^ are different from zero the following equality holds: 

X1AH2 = X^''^ (5.57) 

It is also useful to notice that: 

Ai,4j+2 = -^ + 2f (5.58) 

Again, in analogy with the previous case one introduces the SI2 © Vir modules Q^^-'-* of fixed 
spin j (and tacitly m) built on the highest weights 5.23, and realizes that the x^"^^^ are the 
characters of such modules. 

One then notices that by means of 5.42 and following it is possible to derive the following 
identities: 

Xr,2 = Xl,2r = X^'-'^ (5-59) 
n-1 

X.,2n = 5^x(^-"+''=) (5.60) 

A;=0 

which imply that the modules Vr,2 and Vi_2r are isomorphic, and that in general the following 
decomposition holds (which is in agreement with the results of [6] upon switching to their 
notations): 

n-l 
A:=0 

So that we have finally given a description of the entire Kac table in terms of the modules 
and 

In passing it is very nice to make some simple remarks about W— modules. Notice that the 
W— characters of [29] can be cast in the following form: 



{q)=J2{2j + l)x^'^\q) (5.62) 



X2,i(g)= J2 (2j + 1)X^'^H?) (5.63) 
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implying that the corresponding W— modules have the following structure: 



= © © 2^'"'"^ (5-64) 

^■gl^m=-j 



= © © Q^'"'""^ (5-65) 

j€^N+ m=-j 

and similarly (but remember the different meaning of x''^"'^)- 

XiM= $^(2j + l)x(^^)(g) (5.66) 

X2,2(g)= Yl i^J + m^''\Q) (5.67) 
<2 = © ® 2^''''"^ (5.68) 

^2^2 = © ® Q^'"'"^^ (5.69) 

notice that this agrees with the results of the end of the previous section, since: 

j j 

= n2^''"(Xi) = TT, ® V2,+l,l (5.70) 

m=—j m=—j 

essentially because the sum: 

j 

n^^'^ = 71, ® Ivir (5.71) 

m=-j 

is the projector ttj onto 2j + 1 dimensional SI2 spin j representations times the identity on 
the Virasoro part. 

We observe also that the multiplicity (2j + 1) of each module arises precisely from the mul- 
tiplcity of the allowed values of m for the states |j, m). From these expressions it is also 
transparent that the W— modules are closed under the action of s/(2) rising and lowering 
operators, the action of the diagonal generator can be used instead to twist the monodromy 
of the fermion by a continuous phase, thus generating a flow between the twisted and un- 
twisted sectors. 

A further remark is in order here, the direct sum decomposition for W modules in terms of 
Virasoro modules should be taken literally for (1, s) boundary conditions, however for generic 
(r, s) we shall see from the examples reported at the end of this section that sometimes those 
direct sums should be considered as indecomposable sums. A more formal treatment of these 
cases would then require to introduce the so called projective modules (see for example [18] 
and references therein) in order to describe the case of indecomposable sums. 



34 



5.3 Fermionic form of the BLZ Eigenstates 

We proceed now to describe the exphcit relation between the selection rules and the fermionic 
form of the eigenstates of the BLZ lOM. 

Recalling the lattice selection rules for the vacuum sector, we introduce 2-column diagrams 
of infinite height T) G ^^m' labelled by (Z, r) with both Z, r of length m. 
In general one notices that the following state: 

m 

\V) = \[x^U-X^n^ (5.72) 

i=l 

is such that: 

Lq\V) = w{V)\V) (5.73) 
and one identifies w{T>) as the level of descendance. 

Notice that the states |P) can be brought to a canonical form where modes with the same 
label are coupled by a scalar product, to this goal the following identity proves useful: 

Xm ■ XnXm ' Xl = -\xm ' XmXn ' Xl ,m,n,l <0 (5.74) 

Although the counting of states is correct one has to check that the states |P) are always 
eigenstates of the BLZ lOM. It is possible to check by hand that this is indeed the case up 
to level 6, and it should be true at all levels. 

More generally if we consider a state jP) G Q^^-'^ it will be of the form: 

2j m 

\'^)=U>^-r^..UX-k-X-rP (5.75) 

i=l i=l 

and by using 5.74 together with: 

_1, 

Xm Xn ~^ Xn Xl Xm ~^ X,mXn Xl ^ (5.77) 

it is possible to bring all the expressions to a simple canonical form. 

Indeed, it should be true that for all the modules Vi^j+i, the states jX') are eigenstates of 
the lOM, leaving aside possible mixings due to degeneracy. 

The situation for the modules Q^^-^-* is slightly different, in this case the structure of the 
fermionic state associated to a two column diagram is: 

2j m 
i=l i=l 

The singlet case j = is understood to give rise to the analogue of 5.72. 
In this case the action of Lq is given by: 

Lo\V)= (^w{V)-m-j-^^\V) (5.79) 

This difference is related to the presence of q~"^~^ in the definition of the character 5.56. 
Aside from these differences all the considerations of the previous cases apply also here. 



(5.76) 
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5.4 More on lOM and Symplectic Fermions 

It is well known that Confer mal Field Theories are integrable, since they possess an infinite 
set of independent integrals of motion. Furthermore it is possible [2] to build a family of 
such operators for example by quantizing the integrals of motion of the classical Sine-Gordon 
theory. 

A general expression for these integrals of motion is up to now unknown, but anyway they 
can be obtained constructively as polynomials the of L„ by requiring 

[l2„-i,l2;-i] =0 V/,n = 1,2,... (5.80) 

and that the l2n-i have conformal dimensions: 

{h,h) = {2n- 1,0) (5.81) 

so that l2n-i has spin 2n — 1. 

An expression of the first few of them can be found in [2] [1] and is given by: 

Ii = L„ - ^ (5.82) 



n=l 



T ^T^ ^ + 2r , c(5c + 22) 



12 " 2880 



3 

I5 = ^ ] ^m+n+p, • L^LnLp : +— ^ ^ -^^l-2n-^2n-l + 

|, / n + c c V c + 4 (c + 2)(3c + 20) (5.84) 

c(3c+ 14)(7c + 68) 
290304 

where the : : denotes Conformal Normal Ordering which can be obtained by arranging all 
the Ln in an increasing sequence with respect to n. 

The diagonalization of the lOM has to be carried out inside a Verma Module, and since Ii 
already lifts the degeneracy at different levels of descendance, we are interested in building 
a matrix representation of our l2n-i at a given level of descendance K. 
In general if we choose a level of descendance K and build all the strings of operators 
with Ui > and ^ = we have that, due to the presence of null states, the dimensionality 
of the spanned space is most of the times reduced. 

We have chosen to build the matrix representation of the l2n-i using a sovracomplete set 
of states, and then after diagonalization we identified some linear combinations of Virasoro 
generator as null states, when expressend in their fermionic form. 

The eigenvalues are precisely those obtained from TBA and the relation of quantum number 
to fermionic two column states has been discussed extensively in a previous section. 
In the procedure of diagonalization it proved useful to notice that if | /i + i^) is some J'^-th 
descendant of the highest weight one can obtain a truncated action of the lOM at a 
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given level of descendance, the details of this truncation are reported in the appendix. 
This truncation allows for a safe algebraic computation of the matrix representation of I3, 15 
simply using the commutation rules of the Virasoro Algebra. 

We want now to give a description of how the BLZ lOM can be reconstructed in terms of the 
Symplectic Fermion. In [1] the lOM are defined as the modes of weight zero of appropriately 
regularized polynomials in the stress energy tensor and its derivatives. The first non trivial 
one being 

: T'^iz) : (5.85) 

where normal ordering above (which is not fermionic Wick normal ordering and neither 
simple Virasoro normal ordering for the L^) means essentially something that is defined up 
do addition of the correct multiple of d'^T{z) to the regular part of the T{z)T{w) OPE. 
It may be objected that the comparison with the lOM of [1] is somewhat indirect, however it 
turns out that this comparison can be turned into an exact match by virtue of our knowledge 
of the full spectrum of the lOM in their well known Virasoro form. 

Let us now consider T{z) which is quadratic in the Fermionic fields, one can compute the 
TT OPE by using Wick's theorem for free fields: 

^^^)^^-) - - (7^ + (1^ + (!^ + - - 1^^^^^ - - + 



+ ^{d'T + d'T, + ^-d^T, - Ts){z - wf + 



which is obtained from the propagator: 



(5.86) 



{X''{z)x'{w)) = (5.87) 

It is indeed because of the Wick contractions that one ends up having always a result which 
is quadratic in the Fermionic fields. 
The operators are defeined as: 

T2n =: d--'x{z) ■ d'^-'xiz) : (5.88) 

T{z) = ^T^iz) (5.89) 

where :: now refers to fermi mode (Wick) normal ordering. Notice however that these 
are not the same as [1], which used a different regularization procedure. 
We notice that the operators T2„ together with their derivatives, span certain linear sub- 
spaces in the space of all descendands of the identity operator. The idea is that the BLZ 
lOM lie within special subspaces of fixed weight, their decomposition in these linear spaces 
being fixed by our knowledge of their spectrum. 

It turns out that we can give a decomposition of the lOM which holds at least for (l,s) 
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boundary conditions. 

Expanding the fields in terms of modes one lias: 



where 

and introducing 
satisfying commutations: 



fc-2 fc-2 

Pti = U^rn + 1 + + 1 + j) {5.91] 

j=0 j=0 



N„ = " ^" (5.92) 
n 



[N„, X„] = (5„+m + 5„_m,)Xm (5-93) 

one has that the modes of weight zero of the satisfy: 

S — 1 CO 

{d^T,,{z))^ = l[{2k + j){2Y,nP%N.,) (5.94) 

j=0 n=l 

notice that 

k-2 

nP[% = nl[{{l+jf-n') (5.95) 

j=0 

is a completely odd polynomial in n. 

The operator is indeed a number operator and is action is define on tableaux states by 
counting the number of times #(j) the label j appears in the partitions associated to the 
tableaux: 

N,|l^> = #(j)|^^> (5.96) 

where 

#(j) = 0,l,2 (5.97) 

we now consider operators of the form: 

k 

hk-i - /2T1I = E^"('^'^'~"^^2n)^ (5.98) 

n=l 

now 

2fc-2n-l 00 

(92"T2.(z))^= n i2n + j){2j2^P[t,m^m) (5.99) 

j=0 771=1 

and substituting 

00 k 2fc-2n-l 

hk-i-P2k-^ = Y.mNr,Y.^P^tmCn n (2ri + j) (5.100) 

m=l n=l i=0 
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we further seek to impose for some c„ the condition: 

k 2fc-2n-l 

^2c„(P!':i,^ n (2n + j)) =«fc^''"' (5.101) 

n=l j=0 

for k = 1, 2, 3, . . . one has: 

2ci = ai (5.102) 

2(6ci + C2) - 2c2m^ = asm^ (5.103) 

2(120ci + 20c2 + 4c3) - 2(20c2 + 5c3)m^ + 2c3m^ = agm^ (5.104) 

2(5040ci+840c2+168c3+36c4)-2(840c2+210c3+49c4)m2+2(42c3+14c4)m^-2c4m^ = a^m^ 

(5.105) 

so that one has by solving that: 

00 

Ii - /ri = Y (t2)^ = «i 5Z (5.106) 

m=l 

00 

I3 - /3^'^n = g (92T2 - 6X4)^ = "2 J] m^N^ (5.107) 

m=l 

00 

I5 - = - SOd'T, + 120r6)^ = a.Y^m'Nm (5.108) 

m=l 

00 

(^9^T2 - 126a^T4 + 16809^X6 - 5040X8^^ = "4 J]] m^N„ (5.109) 



17 - /r'l 



m=l 



10080 
and in generah 

00 

771=1 

being 

ak = 2^-''k (5.111) 

notice that the l2fc-i are indeed the lOM, and that their action is straightforwardly diago- 
nalized being expressed in terms of the number operators and simplified through 5.96: 

hk^i\V) = {akY.f'-' + irk'L,)\V) (5.112) 

the fact that the l2fc-i are indeed the lOM is the consequence of the fact that its spectrum 
agrees with the lOM on a complete set of fermionic states and therefore by the spectral 
theorem for hermitean operators they have to be the same object. To be completely honest 
this is simple only in the sectors where the lOM are diagonal, for generic boundary conditions 
one needs to define an indecomposable action of the lOM on couples of tableaux states which 
are logarithmic partners: 

l2fe-i|^^> = l2k^iiV)\V) + \V') (5.113) 
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hk-i\V') = l2k-i{V')\V') (5.114) 

where if (1, r) = V one has that (1', r') = V is such that l'^,r^ G (1, r). In other words V is 
another tableaux allowed by selection rules, made with the same numbers as the tableaux 
V. 

We shall understand from the examples in the next section that such couples of states can 
indeed exist and that the diagonal action of the number operator is obtained when one con- 
siders the logarithmic partner I'D') to be a null vector and that I'D),]!)') belong to different 
fermionic modules Q'^K 

However in this case of indecomposable representations the description in terms of Symplec- 
tic Fermions should be treated with more care, which is not our goal in this place. 

5.5 Examples 

In this section we want to give a comparative description of some Verma modules corre- 
ponding to the same conformal weight. The method we shall employ is direct calculation of 
the matrix form of the lOM at a given level of descendance in the standard lexicographically 
ordered Virasoro basis, we will then compute the Jordan canonical form of such a matrix to 
discover that in many cases it exhibits Jordan Blocks. 

In the cases of (l,s) modules it is well known from the lattice theory that the action of 
the transfer matrix is completely diagonalizable, we shall confirm this observation for the 
quantum transfer matrix D which we can define inpiring ourselves to 3.27: 

D(x) = e^(^)+^^-i c/ni.„-.e(--i)^ 

For a suitable F(x) = F{x)l which is introduced to resemble the structure of the expansion 

3.76, and to take into account 4.2. 

Notice that D satisfies an inversion identity: 

D(x - z^)D(x + z|) = e^(-+^f)+^(--^f )l (5.116) 

In some cases, however, with r > 1 we shall find that the higher lOM exhibit a nontriv- 
ial Jordan structure. Lq, by the way, is always diagonalizable because we are considering 
modules built on the vacuum Q by using strings of fermionic operators. What makes D not 
diagonalizable is the effect of the higher lOM. Again, this is in perfect agreement with the 
lattice theory for which according to 4.7 the eigenvalues of the lattice integrals of motion are 
given by a series of the continuum integrals of motion. The reason why the Hamiltonian on 
the lattice is not diagonalizable in some cases is that it is a superposition of continuum lOM, 
so that even of Lq is diagonal, the lattice Hamiltonian receives contributions from operators 
which are not diagonalizable. 

5.5.1 Vi,5 VS V2,7 

We want to give in this section an explanatory study of the module Vi,5, such a module has 
Ai^5 = 1 and is known to have a null vector at level 5. We shall identify such a state as the 



40 



reason preventing the matrix representations of the lOM to be indecomposable. 

The matrix form of the lOM is found to be diagonal up to level 4, at level 5 one finds that 

the Jordan canonical form of J3 can be obtained by a similarity transformation: 

U-1I3U = J3 (5.117) 

which is explicitly realized by: 
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I3 is written in the standard basis: 

[l_,\1),L_,L_,\1),L_3L_2L_,\1),L_,L\\1),L^_,L_,\1),L_,L^^^^^ (5.119) 
The Jordan decomposition is found to be: 
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(5.120) 



Notice that the size of the matrix is -P(5) = 7 but the dimensionality is known from the 
character to be 6. P{N) is the number of partitions of into as a sum of positive integers. 
It may seem that a Jordan indecomposable structure is emerging for /3 in the module Vi^s 
at level 5, thus one introduces the generalized eigenvectors. 



Pi,t = 0,...,6 (5.121) 

Their virasoro form is simply found by applying the similarity transformation to the vectors 
(1,0, 0,0, 0,0,0), (0,1, 0,0, 0,0,0),.... 

After finding their virasoro form one can go over to the fermi modes, and one has: 

4319 

l3P2 = J^P2+Pl (5.122) 

where 

Pi = 10(72L_5-216L„4L_i-144L_3L_2L_i+108L_3L2i+288L^2^_i-180L_2L^i+18Lii)|l) 

(5.123) 
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is found to be a null vector when expressed in its fermionic form, whereas: 

192 

Po = —X-iX-2 ■ X^3^ + 24xl2X-3 ■ X^i^ (5.124) 
5 

P2 = 27QxUx-2 ■ + 510xi2X-3 • X-i^ (5.125) 

P3 = -25xUx-i ■ X-i^ (5.126) 

P6 = 2700x^6^ (5.127) 

these are the eigenstates to be found inside Q^^\ the other 2 states have to be looked for 
inside Q^^^ 

P4 = 18x^4X^2^ (5.128) 

P5 = 90x^5X^1^ (5.129) 

The discussion of this case is sufficient to show that whenever two states, are degenerate for 
all the lOM they are allowed to mix, although one can always to pick a basis within their 
common eigenspace for which, in this case: 

Po = X-iX^2 ■ (5.130) 

P2 = xUX-3 ■ X-i^ (5.131) 
The case of V2J is radically different. First of all one has: 

V2,7 = Q^'^ © Q^'^ © Q^') © Q(^) (5.132) 

and furthermore the submodule one wants to mod out starts at level 14, therefore the 
dimensionality at level 5 of V2J is precisely -P(5) = 7, for this simple argument pi cannot be 
a null vector. The only candidate for a nonzero pi can be taken from the module Q^^\ so 
that it is natural to suggest (upon suitably normalizing everything): 

Pi = xUxUxt.n (5.133) 

So that in this case the action of I3 (and likewise all the higher lOM) becomes indecomposable 
at level 5. Such a Jordan cell will successive levels of descendance. 

Explicitly at level 6, under the action of L_i one has: 

7]i = L_ipi = 3x^4X^2X^1^^ (5.134) 
which spans a Jordan cell together with some suitable linear combination: 

V2 = aiX^iX-4 ■ X-2 + a2X^2X-i ■ X-i (5.135) 

such that: 



l3V2 = ^^V2 + Vi (5.136) 



at level 7, there are 2 Jordan cells. 
The first one is spanned by: 



ei = (L_2 + 5L\)pi = 21xl,xUxti^ (5.137) 
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and again, some linear combination: 



6 = biX-iX-5 ■ + b2X^2X-5 ■ X-1 (5.138) 

such that: 

/3 6 = ^6 + ^1 (5.139) 
whereas the second Jordan cell is spanned by: 

«i = (-4L„2 + L\)p, = Uxt,xUxt,n (5.140) 
and the usual linear combination: 

"2 = cixlix_4 ■ X-3 + C2xUx~4 ■ X-1 (5.141) 

such that: 

11039 

"2 = -^^Q ^2 + «! (5.142) 

The Jordan cells generated by pi will initially be counted by P{N — 5), but for large 
enough this will change due to the appearence of a rank 3 Jordan cell at level 14. In general 
one will notice that at levels n{2n + 3) = 5, 14, 27, 44, ... a new null vetcor will appear and 
in correpondence one will observe higher and higher rank Jordan cells appearing. One then 
will use the states available from the modules Q^"^^^ to fill up the null vectors spanning the 
Jordan blocks. 



5.5.2 Vi,2 vs V2,4 

Let us now consider a module built on a primary field of dimension A = — |, by standard 
calculations with the Virasoro algebra one can start to work out the explicit form of the 
matrix representation of I3. Already at level 2 things start to be interesting, one finds the 
following Jordan decomposition: 

/ 3367 1 \ 

7 3367 (5-143) 
\ ^ 960 / 

which is obtained by the similarity trasformation: 

U = ( 7 ^ ) (5.144) 

Indeed from the character Xi,2 of Vi^2 = Q^'^'^ we know that at level 2 we have a null vector, 
which corresponds to the combination: 

Po = i-^L^2 + L\)fi = (5.145) 

and a generalized eigenvector which is just an eigenvector because of the equation above: 

Pi = ^^-1/^ = ■ (5.146) 

D 2 2 
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satisfying: 



I3P1 



3367 



-Pi + Po 



3367 
"960 



■Pi 



(5.147) 



960 

The difference witfi the module V2,4 = Q^^^ © QP''^ starts early, since from the character we 
see that it has no null vector at level 2 and since the module Q^'^^ has not enough states the 
only candidate is the highest weight of Q^^-*: 



Po 



Xl3X_lP 

2 2 



(5.148) 



so that we realize that the action of I3 starts to be indecomposable already at level 2. 
The state po will generate the whole module QP'\ which will always appear inside Jordan 
cells having the same multiplicity at a given level of descendance as the corresponding 
module. At high levels however, some Jordan cells will still contain null vectors due to the 
fact that the next module in the sequence which is Q'-^-* is not available for filling up those 
null vectors. The next null vector will be at level 8, and in general one will have that a new 
null vector will appear at level 2n^ for n G N. 

At level 8 the Jordan decomposition of I3 is indeed very big, by the way the interesting part 
is that instead of having P(8 — 2) = 11 rank 2 Jordan cells we have 8 such cells plus one 
rank 3 Jordan cell. This Jordan cell clearly contains the new null vector. 
We reproduce the 6x6 block of I3 which contains such a cell together with a rank 2 cell 
having the same eigenvalues: 
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(5.149) 



notice that the repartition into 1 rank 3, one rank 2 and 1 spare eigenstate sums up to 6, 
precisely as the allowed states from the modules QS^\ QS'^\ Q^^-*. The states spanning the 
cell are: 

(3,1|4,2),(2,1|4,3) e (5.150) 
(3|4,2,l),(2|4,3,l),(l|4,3,2)e (5^15^) 

(|4,3,2,1) G Q(') (5.152) 

At level 10, as it is natural to expect, there are P(10 — 8) = 2 rank 3 Jordan Blocks 
with the same structure as the one appearing for the first time at level 8. It is natural 
to conjecture that each time a new null vector will appear it will bring along a Jordan cell 
of higher rank, so that at level 18 a rank 4 Jordan cell is expected to appear for the first time. 



6 Discussion 

In this paper we have seen how the natural combinatorics describing Critical Dense Polymers 
is the same describing Symplectic Fermions in the continuum limit. Furthermore we have 
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been able to obtain the eigenvalues of the BLZ local involutive charges from Thermodynamic 
Bethe Ansatz. These involutive operators and their eigenstates are described extensively in 
the continuum theory in terms of Symplectic Fermions and the mechanism providing the 
Jordan indecomposable structure of the continuum transfer matrix is discussed for some 
boundary conditions with r 7^ 1. 

A fermionic decomposition of the lOM in terms of the symplectic fermion however is still to 
be achieved in presence of indecomposable representations. Also, the analisys of indecom- 
posable structures for r 7^ 1 boundary conditions should be carried out on the lattice and 
checked against the results presented here in the continuum, in some specific examples. 
This point of view on the integrability of the model allows us to recognize the eigenvalues of 
the BLZ charges in a number of exact expansions holding directly on the lattice, this enables 
us to introduce the involutive charges on the lattice and to identify both their eigenvalues 
and their decomposition on the Temperley Lieb algebra, providing beautiful exact formulae 
for the decomposition of the transfer matrix in terms of conserved quantities. 
This analisys can be carried out for the critical A3 model as well, in this case however it 
would be more natural to decompose the lattice involutive charges on the Clifford algebra of 
7 matrices, the number of generators being related to the system size N. 
Some extension of this work in the direction of arbitrary loop fugacity (3 should also be 
possible. 

Finally we report that shortly after submission it has been possible to extend the analisys 
of this work to the Baxter Q operator as a tangle defined in the Temperley Lieb Algebra, 
as a consequence also nonlocal charges are introduced, and their expressions are discussed 
in [24], which should be read with the spirit of getting some rest after the heavy load of 
50 pages of this publication, and to recover along the way the known local charges on the 
lattice, however with a more powerful and quick method. 
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8 Appendix 

8.1 Bernoulli Numbers 

The Bernoulli numbers B„ are defined as: 




(8.1) 



n=0 



they satisfy: 



5. 



'2n-l — 



,n = 2,3, ... 



(8.2) 
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They appear in the Euler Maclauring summation formula: 

E /(") ~ fmd^ + M±M + f; - r-"(a)) (8.3) 



these numbers satisfy a wide variety of identities, for example: 



Erri^-O (8.4) 



n=0 

can be used for proving identities like: 



2n-l ... B 



2(2n + l)! ' (2n-s + l)!s! 



=0 
2n-l 

4"(2r2+l)! ^ 2^(2n-s + l)!s! ^'^ 

which ensure that the even derivatives of F{t) drop out of the Euler Maclaurin calculation 
3.55. 

Or they appear in the sum of powers used in 3.44: 

fc=l k=l ' 

They are also necessary to go over from the zeta functions appearing in ?? to the expressions 
of the highest weight BLZ lOM 3.40: 

o2n— l„2n 

C(2n) = (-l)""'^^52n (8.8) 
8.2 Proof of an Integral 

We want to give an explanation of how 3.29 was obtained. The integrand is not bounded 
along the imaginary axis and has also double poles. Nonetheless we want to find a way to 
evaluate the integral by means of the residue formula. For this reason we split the double 
poles by inserting a regulator which we also use to introduce a dumping factor along the 
imaginary axis. We will evaluate this integral by using a contour running along the real axis 
and enclosing the poles in upper half plane with a semicircle of infinite radius. 
After introducing the regulator one gets: 

/„(e) = / dx- ^ — (8.9) 



oo 



1 + Y + cosh X 



the roots of the denominator are: 



e2 /. e2,2 



x,{e) = log(-l--±'^(^l + -) -1)~ Z7r(2/ + 1) T (8.10) 
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one then gets by expanding in powers of e and isolating the residues (the divergent part 
drops by virtue of the ± signs) : 

oo 

/„(e) = -Am ^(-z)(-l) V"-i2n(2/ + l)^"-! + 0(e) (8.11) 
/=o 

one then uses in sequence: 



E(^-2^LU.)-LU-.)) 

u,{-i) = {2'~^ - my) 

C(l - 2n) = (-l)3"2^-2"r(2n 



C(2n) 

2n 



to obtain 
so that: 



4(e) = 8^(2^-^" - l)r(2n)C(2n) + 0(e) 



™2n 

4(0) = / dx — = 8n(2i~2n _ i)r(2n)C(2n) 

' 1 + cosh X 



8.12) 

8.13) 
8.14) 

8.15) 

8.16) 



The other similar integral is evaluated precisely with the same techniques. 



9 A Resume on Characters Selection Rules, and Con- 
nections to other Notations 

We define a two column diagram as: 

P= (/i,/2,...,|ri,r2,...) = {l,r) (9.1) 
to each diagram is associated a weight given by 

w{V) = J2h + J2^, (9.2) 

i 3 

then one introduces the sets of admissible two column configurations of maximum 

height M satisfying the dominance U < ri,i = 1, . . . principle of [5], and their associated 
characters, the generalized g— Narayana numbers of [5]. 

E (9-3) 

{J!n)^ = 0, if<„, = (9.4) 
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the finitized characters are then given by: 

X^^(^) = ^-E LTU^ (9-5) 



N-2j 
2 

JV-2 



m=0 



X^(g) = X^^(g) = g^E UZ2,) (9-6) 



n 9 
m=u 

it happens then that we can express everything in terms of the following characters of 
fermionic W modules of fixed isospin j 

X2j+i,i=t'iq) (9.7) 

Q'' = V2,+i, 1 (9.8) 

where Q'^^ is a fermionic module of W = sl{2) ^ Vir spin j (and tacitly projected onto the 

subspace with fixed J° eigenvalue). 

To be more explicit the modules Q^^ are built as: 

Q'^ = {L^n, ■ ■ ■ L.n,X% ■ ■ ■ xt,^mX--,+m+l ■ ■ ■ X'Jl^} (9-9) 

where in this case Xn ^"^^ modes of the symplectic fermion. 

All othe characters can be decomposed as sums of characters of these modules. 

for example: 

Vi, 4,+i = Q'^"' © Q^' (9.10) 

or in a more complicated case 

V2,+i,3 = Q^'^~'^ © Q^^'^ © (9.11) 

Contact with the notations of [18] is made through: 

xl.,Aq) = xT\q) + xM = C^iiU^M + C^li^[N]{q) (9.12) 

suggesting: 

xM - c^^,2^mq) (9.13) 

with similar identifications of Cf'^\N\{q) for even s. 
another identity that should hold is for example: 

X%^,M = C^]2f-ANM + Ca5^[iV](g) + Ci],%ANM (9.14) 
and a lot of similar identities. 
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9.1 Some Virasoro- Fermi Modes Calculations 

We start by analizing the module Q'-'^^ (being also the vacuum module), for which the most 
generic states up to level 6 give, after some straightforward algebra: 

Q (9.15) 

L.2^ = \x-i-X-i^ (9.16) 
L^^n = X-2-X^i^ (9.17) 
(aL_4 + hL\)n = (a + 6)x_3 ' X-i^ + |x_2 ' ^-2^^ (9.18) 
(aL_5 + bL^3L^2)^ = {a + 6)x_4 " X-i^ + aX-a ' (9.19) 
(aL_6 + 6L_4L_2 + cL^3 + dlt^)^ ={a + b + 2c + 3d)x^^ ■ X-i^ + (a + c)x-4 ■ X~2^+ 

(9.20) 

The first few generic states in the module Q^^^ up to level 5 are: 

Xt,n (9.21) 

L^ixti^ = xU^ (9.22) 
(aL_2 + hL\)xti^ = (a + 2h)xU^ (9.23) 

(aL_3 + 6L_2i:_i + cL\)xUn = (a + 26 + 6c)x+4^ + ^(^ - ci)xUx^i ■ X^i^ (9.24) 
(aL_4 + 6L_3L_i + cL\+dL_2L\ + eLiJx^i^ = (a + 26 + 3c + 6rf + 24e)xl5f^+ 

-2' 

(9.25) 



+ ^(-« + c + 2rf)x^3X-i ■ X-i^ + i(« - ^)x^iX-2 ■ X-2^ 



(aL_5 + 6L_4L_i + cL_3L_2L_i + dL^-iL\ + eL%L_i + /L_2L^i + ^L^i)x^i^^ 
+ (a + 26 + 3c + 6c/ + 8e + 24/ + 120^)xi6^ + ^(-a + 4e + 6/)xt4X-i • 
+ (a - c - 2rf)xliX_2 ■ X-3^ + (6 - c - 2rf + e)xl2X-3 ' X-i^ 



(9.26) 

The first few generic states in the module Q*-^^ up to level 3 are: 

xUxU^ (9.27) 

L.^xUxti^ = 2x+3X-i^ (9.28) 
(aL_2 + hL\)xt2xt^^ = (26 - a)xUxt2^ + (2a + 66)x^4X^ifi (9.29) 
(aL_3 + 6L_2L_i + cL\)xt2xt^n = {-a + 12c)x+4X-2^ + (2a + 66 + 24c)xl5X^i^ (9.30) 
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We give also the first few states in tlie module Q^'^^: 

/i (9.31) 

L-ifi = ^X-i ■ X-ifJ^ (9.32) 
(aL„2 + bL\)n = (a + ^)x„3 ■ (9.33) 

^2 2 

^ 3 Qi h 

{aL_3 + 6L_2L_i + cLii)/i = (a + - + -c)x„| ■ + (2 + 4)X-| " (9-34) 

6 3 3 15 

(aL_4 + bL_sL_i + cL\ + dL_2L\ + eLi J/x = (a + - + -c + -rf + — e)x_i ■ 

z z 4 o 2 

, c (i 3 , ,6 c d e, 

+ + 2 + 4 + 2"^^-t ■ ^-f ^ + ^4 ~ 2 " 4 + ' ^"t^-^ ' ^^^^ 

(9.35) 

Tliese formulas, together with table 1 are sufficient to reconstruct the precise coefficients of 
the Fermionic expressions everywhere in the text. 



9.2 Truncated action of the BLZ lOM 

In this appendix we give useful truncations for the action of the BLZ lOM when acting on 
some descendant state at level K: 



, c + 2 c(5c + 22) 
I3 = 2 ^ + Ll- -^Lo + ^2880 ^^"^^^ 

n=l 

K K n-1 2 K 

I5 =3! ( L^nLoLn + ^^(L„(m+n)-^m-^n + L^nL^mLn+m)) + 2 -^l-2n-^2n-l + 

n=l n=2 m=l n=l 

, #>, ll + c 2 c , c + 4 2^(c + 2)(3c + 20) c(3c + 14)(7c + 68) 

+ 1^^^'' 'A~ ' + 576 ^° 290304 

n=l 

(9.37) 
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Table 1. Fermionic structure of the modules Q^"^ up to w{V) = 6, the ** means that the states are degenerate within the same 
module, and therefore linear combinations fall within the same eigenspace. 



